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A NOTE ON THE MATHEMATICAL BIOPHYSICS OF 
AMEBOID MOVEMENTS 


RALPH BUCHSBAUM, N. RASHEVSKY AND HENRY E. STANTON 
THE UNIVERSITY OF CHICAGO 


Theoretical considerations lead to the expectation of some regular- 
ities in the apparently random changes of shape of an ameba. We 
should expect quasiperiodic fluctuations of the ratio of the perimeter to 
the area of the optical cross-section of an ameba, when plotted against 
time. The quasi-period should be characteristic of the species. Pre- 
liminary experiments appear to bear out that prediction. 


This note describes preliminary results of a research problem, 
planned originally on a larger scale. Due to the present emergency 
and the duty in the Armed Forces of one of the authors (R. B.), the 
work has been temporarily discontinued. However, the preliminary 
results obtained appear of sufficient interest to warrant this note. 

A few years ago one of the authors (Rashevsky, 1939) showed 
mathematically that under the influence of the diffusion drag forces 
a cell may exhibit either periodical or quasiperiodical changes of 
shape, by alternately elongating and contracting. It was suggested 
then, that a similar mechanism may underlie the protrusion and re- 
traction of a pseudopod. In an actual case the presence of a large 
number of disturbing factors may obliterate the quasiregularity of 
the phenomenon, and at first glance it may appear to be devoid of any 
such regularity. In the case of a simple elongation and contraction 
of a cell, the ratio of length to width of the cell is the natural quan- 
tity whose variation with time could be studied. In a moving ameba 
such quantities are sometimes not defined at all. The difficulty may, 
however, be overcome. 

An elongation of a regularly shaped cell is always accompanied 
by an increase of the specific surface. If we consider axially sym- 
metric cells, then an elongation also consists of an increase of the 
ratio of the perimeter of the meridional cross-section to its area. For 
a cylindroid cell of length 27, and width 27, that perimeter is equal 


approximately to 
Dea (Ts tts), (1) 
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while the area of the meridional cross-section is equal to 
SAT Ta (2) 
Since for a constant volume V, of the cylindroid cell 


2n7 17.2? = V,. = const. (3) 
therefore, putting 
7 V,/aa—a, (4) 


we have 
p (ghee 3 a (11/%2) +1 (5) 


8 ar, a?!3 (175/19) 2/* 


Hence if 7,/7. fluctuates, so will also p/s, and if for a constant Vo 
we know either one of the quantities p/s, 7 or 11/72, we know the 
other two. 

The quantity p/s has the advantages of being independent on the 
shape and relatively easily determined. In a cell of regular shape we 
may observe 7,/”, directly. In an irregularly shaped cell we may ob- 
serve p/s and define 7, or 7,/7’2 from equations (5). 

If the apparent lack of regularity in the ameboid movements is 
due to disturbing factors, which affect strongly the details of the 
shape, but do not upset fundamentally the mechanism of elongation 
and contraction, we may expect that the quantity p/s will exhibit 
quasiperiodic fluctuations, even though we cannot define such quan- 
tities as length or width. A true periodicity is not likely, since the 
disturbing factors wili affect also the length of the period. By plot- 
ting p/s against time we should, however, expect a fluctuating curve, 
in which the distances between successive maxima or minima, while 
not equal, scatter around a constant value. That value, which may be 
tentatively identified with the theoretical period, should depend on 
the physicochemical constants of the cell, and be a characteristic of 
the species. 

To test this expectation photomicrographs of several amebas 
(Amoeba proteus, Chaos chaos and a limax type) were taken at 
regular intervals of the order of a minute. The perimeters and areas 
were then measured on the enlarged photographs. The perimeters 
were first measured with a map measure, the areas with a plani- 
meter. Subsequently an apparatus has been devised consisting of a 
combined planimeter and a pantograph. The long arm of the panto- 
graph carries a vertical graduated capillary filled with ink and makes 
an enlarged tracing of the cross-section as the point of the plani- 
meter is moved along the original photograph. The amount of ink 
drawn from the capillary measures the length of the perimeter, while 


R. BUCHSBAUM, N. RASHEVSKY AND H. STANTON 63 


the reading of the planimeter multiplied by an appropriate constant, 
gives the area. 


FIGURE 1 


B/s ARBITRARY UNITS 


280 300 LI) Se 260 See 270 350 250 180 210 
TIME IN SEC. —> 


Figure 1 represents a sample of a curve obtained by plotting p/s 
against time for an ameba of the limax type. The intervals between 
the maxima are in seconds—280, 300, 270, 260, 270, 350, 250, 180, 
210, 300, (the last interval is not shown in the figure) with an aver- 
age of 267 sec., and a standard deviation of 47.6 sec. Another similar 
set of measurements, made at a different time, gave a quasiperiod of 
269 sec., with a standard deviation of 85 sec. 

Similar curves were obtained for Amoeba proteus and Chaos 
chaos. The values of the quasi-periods are correspondingly 173 sec. 
s.d. 86.6 sec.; and 9.41 min. s.d. 3.8 min. 

While no definite conclusions can yet be drawn from those pre- 
liminary studies, they seem to indicate that even the motion of ame- 
bas, almost the classical example of something chaotic and random, 
may be described quantitatively and mathematically. 

This work has been aided in part by a grant from Dr. Wallace O. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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It is shown that osmotic coupling between solutes, leads, under cer- 
tain conditions, to differential equations of penetration which have pe- 
riodic solutions. One such set of conditions is that the permeating so- 
lutes have very similar colligative properties, but that their concen- 
trations in the external medium be quite different. It is noted that iso- 
topes may satisfy these conditions exactly, and that they have been ob- 
served experimentally to penetrate as predicted. The differential equa- 
tions bear a resemblance to the interactive equations of Volterra and 
Lotka. 


Certain reaction-systems have from time to time been proposed 
to account for at least some observed biological periodicities. The 
analysis of one of these has been given quite completely by A. Wein- 
berg (1938, 1939) and N. Rashevsky (19388, p. 54), namely the cou- 
pling 

Ac, /dt = A410, + AyeCe , 

AC2/ dt = M10, + Agele , ~ 
which has the physical interpretation of mutual catalysis without 
combination. A second system was proposed by A. J. Lotka long ago 
(Lotka, 1910, 1920), 

dc,/dt = (dy, + Ay2Ce) , 


deo/dt = Cz (G2 ++ A101). 


Originally, Lotka gave the system (2) the interpretation of two con- 
secutive irreversible reactions in which the first substance catalyzed 
_ its own decomposition. In this guise, Lotka realized that the connec- 
tion with reality was quite remote.* At least to some extent, reaction- 
system (1) is subject to the same criticism as Lotka’s original hy- 
pothesis. Catalytic systems seem invariably to feature the forma- 
tion of intermediates, and their analysis is plagued by corresponding 
difficulties. 


be it., p. 274. Later, however, in the hands of Lotka and Volterra, the 
Pin (a tek on the highly significant interpretation of interaction between 
predator and prey species, and in this form it was experimentally verified by G. 
F. Gause. 


(2) 
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In attempting to study the simplest cases of penetration of mem- 
branes under general thermodynamic conditions, the author was led 
to equations which in a sense are referable to equations (2), and 
which in certain cases may explain periodic or oscillatory penetra- 
tion. These relations seem justifiable by elementary thermodynamic 
reasoning. 

Consider the binary system of substances 1 and 2. For the 
chemical. potentials of these components we have, 


fa = wi it RT log ¢, + arCe ; 

[2 = Ls ae RTlog Co FP OleC1, Be 
where a, and a are constants, and the remaining symbols have their 
usual meanings, 

The inclusion of the terms a2¢, and a,c, may be argued in the fol- 
lowing way. The expression for u, contains a term which represents 
the work done in bringing one mole of substance 1 from infinity into 
the solution. This work in general depends on the presence of the 
other solute, substance 2, and therefore is a function f(c.) of ¢.. For 
not too large values of c. we may expand f(c.) and retain only the 
linear term, a, ¢.. The constant term of the expansion is included in 
u’,. A similar argument applies to a,c, . 

From equations (3) we have for the force of diffusion on, say, 
the first component, 


> 


RT ; 
erN came ores eee gO NAEY (4) 


whence the current of diffusion is, 


4, = D,¢,F,= —D,(RT V 6, + 1 4 V C2). 


and the penetration over the vector surface, s , 1s 


> 8S 


S-i4,=D,¢,F,-S=—D, S..(RE Wieut che ayecae (5) 


The quantity D, is not the diffusion coefficient, but is proportional to 
it. For the sake of simplicity we henceforth consider motion parallel 
to the z-axis only. We then have the single scalar equation, 


ae Q 
Se tg = —P), (RT ile De oe Oy C1 ae Sz) ’ (6) 
Ox Ox 


with S, as the area normal to the z-axis. Let us imagine that the 
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membrane to be penetrated is a homogeneous and isotropic barrier 
of thickness dx and of a phase distinct from the two regions which 
it separates. These two regions, however, (corresponding to exterior 
and interior of the cell) we regard as being the same type of phase, 
to follow the customary procedure. Let the partition coefficients be- 
tween the membrane substance and the two regions be p, and pz re- 
spectively for substances 1 and 2. Under these circumstances we can 
approximate 0c,/dx and 0c./dx across the membrane. (See Figure 1.) 


$x 
CIG= Gc, 
CICa CalC, 
FIGURE 1 


Taking, for instance, ¢°, | ec’, as the concentrations at the immediate 
exterior-membrane interface, and c’, | ¢, as the concentrations at the 
immediate membrane-interior interface, we may write approximately, 


Ci C's — C's 01 — 
ta ee as ng a Pe aT ra 
0x membrane 5x 64 


and on substituting into (6), 


(c%, — G1) (c% — C2) 
(Sz 5) membrane —— == Ds Sz (re, eis ees = Oy Ci Pro eee es ; (7) 
8x 8x 


or 


—D,S,RT9, —O1 C1 Pr DS, 
(Sz tz) memb. = aa (CoC) =F ae pee (cle) a 23) 


It is seen that the doublet, 


—D, S, RT: —, D2 D, S, 
eee }= Gansta) (9) 
8x bx 
rather than the usual h,, , characterizes the permeability of the mem- 


brane. 
Combining equations (4), (8), and (9), we are led to the two dif- 


ferential equations, 
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de, /dt = hy, (6°, — C1) + Ire CC's = 265); 
dc. /dt = ical Co ited c) =f Noy Cy (c% ro C3). 


Assuming that the external concentration of these substances is sen- 
sibly constant, we may write, 


(10) 


Ac, /dt = My, Cy + Ayy Cy Cy + Arg , 


(11) 
AC2/ At = Any Co + Ag, Co Cy + Mos , 
with 
Ohi Mg Cy = Nig ls Nee. 
Aye = —hiz , Oey = <her ’ (12) 


Gis = Ny 6s , og = Ngp Co . 


Let us set uw, = ¢, — Gist, aNd Ue = Cy — Aggb. We may then continue 
formally by substituting into equation (11), and obtain: 


du, /dt = Ay, Cy + Aye Cy Ce , 

(13) 

AUsflt = Age Co + Ge; Cy Co. 

Eliminating t by division and rearranging we have, 
(Qe2/C, + Qn) Ay = (Agr /Co + Aya) dUs « (14) 

But in view of our change of variable, this expression is equivalent to: 


(Ge2/C, + M1) dey + (Ayr /Co + Aye) Aosdt 
= (Aq2/C, + Ay1) Ay3At + (Ay /Co + Aye) ACy , 


whence, integrating formally, we have: 


Mog lO Cy + My1€C, — (Ary LOG Cy + Aye Cy) 


15 
= K + Gy.Q13 is dt/C, — Ay,Ass fe dt/c. + (Gy3Mg1 — AyMy3) TF, oo 

with K as a true constant of integration. The equation 
Gog log C1 + Ae, Cy — (@41 log Cg + Gy. C2) = K (16) 


is the integral of the system (2) arrived at by Lotka and Volterra 
(Lotka, 1910, 1920, 1925, Volterra, 1931, 1935), and shown graphi- 
cally by these authors to represent a closed curve, thus indicating wun- 
damped oscillations in time of c, and c.. It is evident that equation 
(15) can be studied from the point of view of regarding equation 
(16) as an equation in which K is a function of time. 

First, however, it is necessary to consider a simple analytic proof 
showing that equation (16) represents, indeed, a closed curve. Take 
the state of affairs when c, has a definite value, say €. We may write, 


MANUEL F. MORALES 69 


K + dy, log & + Gy. & 


Aes 
¥Y(¢c,) =log ce +—aq= (17) 
Ass Ags 
(c,) has one and only one maximum when, 
sees + Qo1/Ax2 = 0, or CC, = = ao / Oe 3 (18) 


(since for all c, , d?W/dc?, < 0); 


whence if there are any solutions at all to (17) there will be one or 
two, according as, 


oo _K+ a, log é+ ay é 
102 — {| = ——__., 
Ae 


(19) 
Azo 


If the inequality of (19) holds (as it can be made to do, by proper 
choice of &), then the line c, = & removes two points from the set of 
points represented by equation (16), and it follows from the simple 
relations (18), that the remaining sets are connected. There are there- 
fore two partitions, whence the set (16) is homoeomorphic to a circle, 
i.e. itis a Jordan curve. (See, for instance, Newman, 1931, page 54). 
In view of the relation (19) it follows that the roots of (17) 
will be far apart or close together according as K is small or large 
respectively, or more pictorially: as K grows larger, the integral 
curves represented by equation (16) are reduced in their vertical di- 
ameters (and conversely). Let us now return to equation (15), in 
which “K” is, in effect, a function of time. Successive segments 
(traced out in time) of the ¢,—-c, curve represented by equation (15) 
can be regarded as segments of integral curves of (16) whose di- 
ameters are being continuously reduced, provided the expression, 


Az2M43 ab: dt/C, — A123 ike At/Cy + (y3Ae, — AyoM23) t (20) 


increases in time. Under such conditions the curve (15) will have 
the general property of a spiral winding inward, which implies that 
C, and G are both approaching steady values in time, but approaching 
them in oscillatory fashion (damped oscillations). We cannot go fur- 
ther with (20) in general terms, except to remark that the two inte- 
grals will be determinate throughout, for the concentrations pass 
from some non-zero value to another non-zero value (steady state). 
Let us now see if the increase in time of expression (20) can be 
assured in some specific case. To do so we must examine the nature 
of the constants, a;;, in the case to be studied. If substances 1 and 2 
are two radioactive isotopes, or at least two rather similar substances, 
we might reasonably expect that the permeability constants, h , would 
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be of the same order of magnitude. In the isotopic case, it would be 
proper to equate them, since such properties as diffusion velocity de- 
pend only slightly on nuclear structure. Let us take h,;; as the princi- 
pal permeability and hi; as the cross permeability, and let us suppose 
that the integrals are of the same order of magnitude, say I (as would 
be the case if the penetration rates were similar). We may then write 
expression (20) as: 


— h?T (0s — C's) — hishiz (Cs — C2) tb, 


whence if c°, < ¢% the conditions for the increase of K with time are 
clearly satisfied. A most suggestive set of experiments, in this con- 
nection, is that of S. C. Brooks on the uptake of radio-potassium by 
Nitella (Brooks, 1940). Normai potassium is of course penetrating 
simultaneously, and the external concentration of normal potassium 
is distinctly greater than that of radio-potassium. The scrupulous 
care with which such experiments were done argues for the reality of 
the oscillations observed; and the striking similarity between the ob- 
served phenomenon and the rough qualitative predictions given above 
suggests that the experiments might well be explained by this type 
of coupled penetration. A more extended analysis of these equations 
—which it is hoped may be carried out in the future—is necessary be- 
fore application to other than isotopic cases. 
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ON THE THEORY OF CELL DIVISION FROM THE POINT OF 
VIEW OF THE PRINCIPLE OF MAXIMUM 
ENERGY EXCHANGE 


HENRY E. STANTON 
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An analysis is made of the equation of cellular elongation during 
division, which was derived by N. Rashevsky from the principle of maxi- 
mum energy exchange. The method used is the same as that employed 
by H. D. Landahl in discussing a similar equation deduced from the 
theory of diffusion drag forces. The differential equation is expanded in 
a power series of the relative elongation, and in this way is reduced to 
the form studied by H. D. Landahl, which has been shown to agree very 
well with experimental data. An estimate of the order of magnitude of 
the universal constant 7, which appears in the generalized Hamiltonian 
principle, is made, and 7 is found to be of the order of 10-4 sec. 


Some time ago, N. Rashevsky (1948a) suggested a generalized 
- Hamiltonian principle, applicable to biological phenomena. The usual 
Hamiltonian principle of physics is obtained as a particular case of 
the generalized one. As an example of a possible application, Ra- 
shevsky has treated the problem of cell elongation during division. 
The differential equation for cell elongation obtained in this manner 
looks very similar to the one obtained previously on the basis of dif- 
fusion drag forces (Rashevsky, 1940; Landahl, 1942b), and can be 
shown to give integral curves of a similar shape. The purpose of this 
paper is to investigate in detail to what extent the new equation of 
Rashevsky’s agrees with data obtained by R. Buchsbaum and R. R. 
Williamson (1943) for Arbacia eggs, which were correctly predicted 
by H. D. Landahl (1942b) from the old equation. We shall see that 
Rashevsky’s equation can be reduced to the same form to which Lan- 
dahl reduced the old one. 
Denoting by 


7, = one half of the “length” of the cell 


1) = radius of the cell when spherical 


c, = concentration of the single metabolite at great dis- 
tances from the cell 


D = internal diffusion coefficient 
71 


72 ON THE THEORY OF CELL DIVISION 


V = total volume of the cell 

w= heat of reaction per gram of the substance pro- 
duced 

a= constant of proportionality such that the rate of 
production of the substance:is equal to a¢ with 
é the average internal concentration of the me- 
tabolite 

r= universal constant in the Hamiltonian principle 

n = average coellicient of viscosity of the cell 

y = surface tension 

S= total surface of the cell 


Rashevsky’s (1943b) elongation equation for a cell may be written 


dr,  2VD wareor 8 — T°) Ts 

dt un [8D (2r2 + 73/7) arr, |? 
Ar’y dS (1) 
9V 7» dry 


The equation discussed previously by Landahl (1942b) is 


T+ dr, 


1, — To at 


_ RTqur’ 1 yf te \ 
4M yD. 1+ (1 +b)1,3/2r,3 + br,8/2/1,3/2 eS st % ) * (2) 
where 
b=2Di/De3 V=4a7,72/3=4273/3, (3) 


with D; and D, being the internal and external diffusion coefficients, 
respectively. 


Introducing relation (3) into equation (1) and putting 
r=T,, W=ar?/6D, K,=Vweres/18 27 D (4) 

we find 
O pe ee ane oe 
dt re[1—a'ro/r + (ro/r)?/2]? 38a To dr’ 


Similarly, equation (2) may be transformed into 


(5) 
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ar 73/2 — 73/2 


— 


SS Oe ee ee ee 
dt PATIL 0 re[P yee (1b) i/7)*] 


y 3/2 =e TF (6) 
- 27 13/2 
where 
K,.=RT qur?/4M7D,.. (7) 


In equation (5) the denominator in the first term is always posi- 
tive if a’ < 3/2 and consequently, the parameters in equation (4) will 
be restricted according to this condition. The similarity is now ap- 
parent between equation (5) and equation (6), provided dS/dr has 
the proper functional form. 

For the determination of the explicit form of dS/dr , Rashevsky 

(1948b) has suggested two approximations: 


I. Approximate the cell shape by an ellipsoid, 
S=4a7) +16272(7r/% —1)?2; dS/dr=382a(r—1%). (8) 


II. Approximate the cell shape by a cylinder, 
S215 (Tor) 2 5-dS/dr =n27(To/7) 2. (9) 


As a third alternative, one might assume that a dividing cell can be 
represented by two intersecting spheres whose radii are equal and 
whose common interior is counted only once. During the elongation 
process the centers of the spheres will separate and the radii change 
in such a manner that the total volume enclosed by the composite sur- 
face will remain constant. F 

Calculating the volume and surface of the two spherical segments 
constituting the cell, and remembering that V,S,7, 7, refer to the 
whole cell as defined in the expressions above, it follows that 


V/2—a27° (hk — 7/3), S=2 Vif ree Aas 5 


where R is the radius of each spherical segment. Since V = 427,°/3, 
we may solve for S in terms of r and 7. The result is 


III. S=4a(27re+7%)/3r; dS/dr=8a(r? —7%°)/3r?, (10) 


the last expression being obtained by differentiation. This shape is 

a closer approximation to the shape found in the biological phenome- 

non of cell division in the intermediate stages, but the final results of 

the calculations in this paper are not materially changed by it. 
Following Landahl (1942b), put 


r— 1 
: (11) 


e = 
To 
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1f either expressions (8), (9), or (10) and expression (11) is substi- 
tuted into equation (5), and the right side of that equation is expanded 
in a power series in e, the following form of the expansion 1s found: 


de/dt=e(A — Be—Cé+ De®?) +E (12) 


where the coefficients A, B, C, D, and E have somewhat different 
forms depending on which of the approximations I, II, or III are 
used for dS/dr. They are given by the following relations: 


12K, 3.29 Hal 6.49 


I. = —_———__ — ; — ; E=0; 
r(3—20')?  3n% 37 To bias 
AK, (15 + 2a’) 2 3.2 a AK, (21 + 2a’) 
= 7 9(3 — 2a’)? 3% f,(3'=— 20)" 
12K. 
Il. Az oe ee t ; B=y/81%o; 
9 (38.— 2a:)3 21 
(14) 
E=- bs : 
37 To 
__ 4K, (15 + 2a’) ped __ 4K, (21 + 20’) 
(8 —2a’)® 4807” 1) (3 — 2a’) 
ev sand, 
128 y 7 
12K 8 8 
1 fagh erator eh esa B=—_; E=0: 
7) (38 — 20)? 37%, 3 To 
__ 4K, (15 + 20’) % 8 y ee AK, (21 + 2a’) ioe 
-#4(8 — 2a)? 97,” é (3 — 2a’)? * 


As Landahl (1942b) has shown, equation (6) also reduces to the 
form of equation (12), but with FH =0. Not as many terms were re- 
tained in Landahl’s calculations, and since in cell division, the maxi- 
mum value which e may have is 0.6, the term involving e* may be 
dropped except for values of a’ very close to 3/2. 

Equation (12) has been shown by Landahl (1942b) to represent 
very well the data obtained experimentally for the elongation during 
division of Arbacia eggs. Hence, equation (1) derived from the gen- 
eralized Hamiltonian principle is shown to agree with those data to 
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the same extent as equation (2), from which the data were actually 
predicted. 

The expansion (12) together with numerical values of the co- 
efficients A, B, and C as calculated by Landahl (1942b) for Arbacia 
eggs makes it possible to evaluate approximately the universal con- 
stant + which enters into the generalized Hamiltonian principle (Ra- 
shevsky, 1943a). The average values obtained by Landahl are 


A=0.023; B=0.052; 7=8.4X10?cm. (16) 


It is also necessary to assume plausible values for many of the 
other physical and chemical quantities listed at the beginning of this 
discussion. Some of them have been evaluated by Landahl (1942b), 
and others of a more purely physical and chemical nature have been 
determined experimentally. It is first necessary to examine numeri- 
cally the magnitude of the constant a’ as defined in equation (4) asa 
check on the validity of the expansions. Considering that for Arbacia 
eggs the total reaction rate aé is of the order 10-7 gm.cm.- sec. (Ra- 
shevsky, 1940), we have with ¢ ~ 10 gm.cm.?, a © 10° sec. As 
a biologically plausible value for D we shall put D ~ 107 cm.” sec"? 
(Rashevsky, 1940). We then have a ~ 10°, and consequently, the 
expansion is valid, and furthermore, we may completely suppress the 
quantity 2a’ in comparison with 3 in the equations of the coefficients. 

If equation (18), (14), and (15) are solved for the constant 
K, in terms of the coefficients, the results for the three cases are 


I K,=0.75 7 (A + B/2) =12.5 X 10°, 
II K,= 0.75 7,.(A + 4B) =59 X 10°, (17) 
Il K,=0.75 7, (A + B) =19.1 X 10° 


the numerical values following from equation (16). 

Having evaluated the constant K,, we may now calculate the 
value of the universal constant 7 from equation (4), although plau- 
sible values must be assumed for some of the constants entering into 
the relation. The volume V of the cell can be calculated from relation 
(3), since 7, is known, and has the value 1.64 < 10-7 em. The exter- 
nal concentration c is higher, but of the same order of magnitude as 
é, hence, we may take ¢@ ~ 10“ gm.cm.*. The heat of reaction w in 
ergs per gram may be estimated from known heats of reaction for a 
large variety of organic compounds. For the pure hydro-carbons, w 
does not change widely for compounds with large numbers of carbon 
and hydrogen atoms (© 50 atoms to the molecule). If the molecule 
includes oxygen atoms, the heat of reaction is generally lowered some- 
what depending on the relative proportion of oxygen in the molecule, 
and consequently, an average value must be assumed for the quantity 
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w. We shall take w ~ 10" ergs per gm. 
Solving equation (4) for 7 we obtain 
18”2DK 
Parte enthrietremse ne irk) Gi (19) 
V wa? Cy 


and for the three approximations, 


ip 7 = 4.3 X 10% 7 sec. 
II. 7= 20.3 X 10+ 7 sec. (20) 
III. 7= 6.6 X 10 7 sec. 


Following Landahl’s discussion (Landahl, 1942c), we may take 
as a plausible value 7 ~ 107+ cm. sec.. This gives for z the order of 
magnitude 

7 © 10+ sec. (21) 


In view of the uncertainty of the values of the constants used in 
evaluating +, the value above must be taken with great caution. If 
7 is of that order of magnitude, then in order for the “energy flow” 
term in the Hamiltonian to become comparable with the regular 
“energy” term, the changes in energy flow must be at least a thou- 
sand times as large as the changes in energy, when the configuration 
of the system is altered. In the inorganic world such systems have as 
yet not been observed. 

The author is greatly indebted to Dr. N. Rashevsky for suggest- 
ing the problem and for many helpful comments. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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A NOTE ON NEURAL NETS 


H. D. LANDAHL 
THE UNIVERSITY OF CHICAGO 


_A neural mechanism is described which provides for a perception 
which is invariant with respect to movement of the stimulus. 


Consider a net in which a neuron N is stimulated equally by each 
of a large group.of equivalent neurons whose afferents are spread out 
uniformily over some surface. Let some part A of this area be sub- 
jected to a suddenly applied constant stimulation. If the surface den- 
sity of the afferent endings is 6, then the total stimulation o(t) act- 
ing upon N at time ¢ measured from the application of the stimulus 
will be given by 


a(t) —A o[z,€% 19,6" + o(1— 6%) — p(1 = e**)] (1) 


where ¢ and wy are constants, provided that all the afferents have the 
same values of « and 7 initially. 

Suppose now, that instead of the stimulus being applied to a fixed 
group of the afferent neurons, the stimulus is moved over the area in 
any manner subject to the condition that the total number of neurons 
stimulated is constant. This requirement is essentially the same as 
requiring only that the area stimulated remain constant. Any move- 
ment of this sort can be broken down into a process in which a neu- 
ron which has been stimulated is removed from the stimulated group 
just as an unstimulated neuron is added to the group of stimulated 
neurons. Then we will have shown that o(¢) is independent of any 
movement of the stimulus if it can be shown that when two equivalent 
neurons are stimulated in immediate succession, the total « produced 
at any instant by them is the same as that due to either one being 
stimulated continuously. 

Thus, consider two neurons, the first stimulated in the interval 
0 <t <7 , and the second in the interval 7 < t, 7 being arbitrary. If 
each neuron has initial values «, and j7,, then for0 <t< 7, 


th) ene ee err Ole Cy op (Le); (2) 
while fort > 7,. 
Tt 
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o(t) =o(r) e289 bp 1 68) ee (3) 


But equations (3) and (2) are identical, as can be shown by simple 
substitution of o(r) from (2) into (3), and thus the argument is 
completed. 

If the stimulus, instead of being constant varies with time, one 
can obtain an expression for o(t) due to the two neurons, in the in- 
terval t > 7, the expression (2) plus the term 


aes ft e*[b(t— 7) — d(t) ]dé 


and less an analogous term with a and ¢ replaced by b and y. These 
terms vanish only if the stimulus is constant or periodic and of pe- 
riod +. Thus, if a number of neurons are arranged in a line, then 
movement would produce no change in a(t) only if the stimulus is 
constant, and a constant movement would produce no change in o(t) 
if the stimulus is periodic with time and moves with the particular 
velocity given by the interneural distance divided by the period. 

The bearing of the above on the theory of the “Gestalt” is ap- 
parent. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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EQUIVALENCE OF THE THEORIES OF NERVOUS 
EXCITATION OF HILL, MONNIER, 
AND RASHEVSKY 


ALSTON S. HOUSEHOLDER 
THE UNIVERSITY OF CHICAGO 


New and direct proofs of the equivalence of the three theories of 
nervous excitation are developed, and by the introduction of the Haupt- 
nutzzeit as the unit of time and a special unit for the measurement of 
the excitatory state, a normal form for Rashevsky’s equations is obtained 
which exhibits clearly the relations among the various independent para- 
meters. 


The equivalence of Hill’s and Rashevsky’s theories of nervous ex- 
citation has been demonstrated in several ways by Offner (1937), 
Young (1987), and Householder (1939). The following demonstra- 
tion seems to possess certain advantages of simplicity and directness. 

If we write x, and x. (or y, and y.) instead of ¢ and 7, the most 
general two-factor linear theory of excitation may be written, in ma- 
trix form, 

to— OAL — Soyer OS, (1) 


where x’ = dx/dt , x , x° and a are column vectors, ais a2 X 2 matrix, 
and S, the measure of the stimulus, is a scalar. Excitation occurs 


when 


Taba (2) 
and if 
r—cy, (3) 
the transformed equations 
y=biy-y)+fS (b=cac,b=c'a) (4) 
represent an equivalent theory provided the relations (2) and 
Y1 = Ye (5) 
are simultaneously satisfied. Evidently, c must have the form 
i ay SOL Ns Ke 
o( r Bea teak (6) 


as is evident when equation (3) is written in scalar form. But if the 
79 
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matrix b is to be a diagonal matrix (Rashevsky’s form), then the 
diagonal elements of b must be the characteristic roots — 4, and — fe 
(necessarily negative) of the matrix a, and the columns 7” and r? of 
ce must satisfy 


at =— mr, (7) 


as is well-known and easily verified. We may take 7 to be any posi- 
tive number, after which choice the transformation is unique. 

As a possibility for the further normalization of the equations, 
we note that they now have the form 


Yi tu(yi- ye’) =BS. (8) 


A change in the unit of time 
= ah b 


replaces yw; and f; by Tu; and T$;. A natural unit of time is that 
required for rheobasic stimulation. If this is assumed in equation 
(8), it follows that, since stimulation occurs at t = 1 (in the units 
of time now employed), 
Bi / Bo = eh fer 
or 
pi=xemt, Bp. eb, 


We may then choose the units for y so that x = 1 and have 
Y's + mi(yi — Yd) = Ser. (9) 


Since the measure of the y enters only through the ratio of the differ- 
ence Yy.° — y,° to a measure of some particular form of stimulation, 
say the rheobasic constant current, the four intrinsic parameters are 
exhibited as this ratio, the unit of time, and the parameters u;. Hill 
introduces the further restriction of assuming normal accommodation 
which reduces the pardmeters to only three and we are able to write 


A= log, B= log dA, 


; 10 
File b/ Me. Bae 

The reduction of Rashevsky’s equations to a single second-order 
equation (Monnier’s form) is most easily effected if we write « and 
7 in Place of Rashevsky’s « — e) and j — j,. Rashevsky’s differential - 
equation may be written in operational form 


(D+k)e=KI, 


(D=d/dt). (11) 
(D+m)j=MI, 
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With a suitable choice of the units for e and 7, the condition for ex- 
citation may be written 

o— 9 el (12) 
If we operate upon the first of equations (11) by (D + m), upon the 
second by (D + k), and subtract, we have 

Cine Z 5) 


d?o/dt?+ (k+m)do/dtt+kmo 
=(K—M)dI/dt+(mK-kM)I. 


To obtain the initial conditions, o) and a)’, we note that a discon- 
tinuity in J(t) in amount A I results in proportional discontinuities 
ine’,j',ando but not ine,7,oro. The discontinuity in o’ is easily 
found to be 


(13) 


Ac =(K—M)AI. (14) 


This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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JOSEPH NEEDHAM. Biochemistry and Morphogenesis. 1942. xvi + 787 pp. Cam- 
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The vegetative activity of cellular aggregates has received but scant con- 
sideration from the point of view of mathematical biophysics, Rashevsky’s theory 
of gastrulation being the only instance. The theory of organic form, as developed 
by the same author, contains no reference to the cellular constitution of the or- 
ganism, and provides no mechanism for the attainment of the form under con- 
sideration, but deals with the end result only. It is to be hoped that some theo- 
retical study may be made in the near future of the physical mechanism involved 
in morphogenesis. \ 

Needham’s book, Biochemistry and Morphegenesis, should prove very valu- 
able in providing a stimulus and a background for such studies. While it is true 
that the author’s interest attaches chiefly to the chemical processes, and while 
there is as yet little apparent physical significance to be found in some of the 
problems treated, nevertheless, with even a casual reading of the book one finds 
many places where quantitative theoretical formulations seem called for, and 
ultimately, chemistry and physics are bound to coalesce, in morphology as else- 
where. From this point of view, the molecular and electronic structures of the 
various inductors will take on increasing significance as progress is made in the 
field of morphology and also in the field of theoretical chemistry. All this the 
author strongly suggests at many points, especially in his discussion of polarity. 

The central theme of the book concerns the manner in which the morpho- 
genetic inductors act upon the reacting elements of the developing organism. 
Thus the book is divided into three parts: The Morphogenetic Substratum, The 
Morphogenetic Stimuli, and The Morphogenetic Mechanisms. The first deals 
principally with the chemical constituents of the egg, the nutrient material, and 
the maternal organs involved in generation and development. The greater part 
of the third deals with the biochemistry of respiration and metabolism, but opens 
with a discussion of the dissociability of the processes of growth, differentia- 
tion, and metabolism, includes a section on heterauxesis and closes with a dis- 
cussion of polarity. Heterauxesis, incidentally, is the term which the author 
recommends for referring to the relations between the rates of growth of dif- 
ferent parts of the same organism, relations which are representable by a recti- 
linear trajectory on a double-logarithmic plot in a very great many instances. 

While the third part represents, in a sense, the culmination of the book, the 
main body of it is contained in the second part. This gives an extensive treat- 
ment of the experimental work dealing with morphogenesis, as well as a critical 
evaluation of the terminology employed by various workers. This second feature 
is hardly less valuable than the first, since a carefully constructed vocabulary 
contains already the framework of the science. And if such terms as “compe- 
tence”, “determination”, “potency”, and the like, are so clearly in need of quan- 
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titative, rather than qualitative, treatment that a feeling of impatience is hard 
to suppress, it is a virtue rather than a defect in the book that the need for quan- 
titative developments should be thus clearly—if not explicitly—indicated. 

It is not necessary to comment on the important role that must be played 
in the development of the subject by such a comprehensive outline by an out- 
standing authority in this field. Mathematical biophysicists in particular, will 
find it invaluable when they do undertake any work on problems of development. 


ALSTON S. HOUSEHOLDER 
The University of Chicago 
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